Demographic inferences and climatic niche modelling shed light on the evolutionary history of the emblematic cold-adapted Apollo butterfly at regional scale.

Cold-adapted species escape climate warming by latitudinal and/or altitudinal range shifts, and currently occur in Southern Europe in isolated mountain ranges within "sky islands". Here, we studied the genetic structure of the Apollo butterfly in five such sky islands (above 1,000 m) in France, and infer its demographic history since the last interglacial, using single nucleotide polymorphisms (ddRADseq SNPs). The Auvergne and Alps populations show strong genetic differentiation but not alpine massifs, although separated by deep valleys. Combining three complementary demographic inference methods and species distribution models (SDMs) we show that the LIG period was highly unfavourable for Apollo that probably survived in small population in the highest summits of Auvergne. The population shifted downslope and expanded eastward between LIG and LGM throughout the large climatically suitable Rhône valley between the glaciated summits of Auvergne and Alps. The Auvergne and Alps populations started diverging before the LGM but remained largely connected till the mid-Holocene. Population decline in Auvergne was more gradual but started before (~7 kya vs. 800 ya), and was much stronger with current population size ten times lower than in the Alps. In the Alps, the low genetic structure and limited evidence for isolation by distance suggest a nonequilibrium metapopulation functioning. The core Apollo population experienced cycles of contraction-expansion with climate fluctuations with largely interconnected populations overtime according to a "metapopulation-pulsar" functioning. This study demonstrates the power of combining demographic inferences and SDMs to determine past and future evolutionary trajectories of an endangered species at a regional scale.

Altitudinal Zonation of Green Algae Biodiversity in the French Alps.

Mountain environments are marked by an altitudinal zonation of habitat types. They are home to a multitude of terrestrial green algae, who have to cope with abiotic conditions specific to high elevation, e.g., high UV irradiance, alternating desiccation, rain and snow precipitations, extreme diurnal variations in temperature and chronic scarceness of nutrients. Even though photosynthetic green algae are primary producers colonizing open areas and potential markers of climate change, their overall biodiversity in the Alps has been poorly studied so far, in particular in soil, where algae have been shown to be key components of microbial communities. Here, we investigated whether the spatial distribution of green algae followed the altitudinal zonation of the Alps, based on the assumption that algae settle in their preferred habitats under the pressure of parameters correlated with elevation. We did so by focusing on selected representative elevational gradients at distant locations in the French Alps, where soil samples were collected at different depths. Soil was considered as either a potential natural habitat or temporary reservoir of algae. We showed that algal DNA represented a relatively low proportion of the overall eukaryotic diversity as measured by a universal Eukaryote marker. We designed two novel green algae metabarcoding markers to amplify the Chlorophyta phylum and its Chlorophyceae class, respectively. Using our newly developed markers, we showed that elevation was a strong correlate of species and genus level distribution. Altitudinal zonation was thus determined for about fifty species, with proposed accessions in reference databases. In particular, Planophila laetevirens and Bracteococcus ruber related species as well as the snow alga Sanguina genus were only found in soil starting at 2,000 m above sea level. Analysis of environmental and bioclimatic factors highlighted the importance of pH and nitrogen/carbon ratios in the vertical distribution in soil. Capacity to grow heterotrophically may determine the Trebouxiophyceae over Chlorophyceae ratio. The intensity of freezing events (freezing degree days), proved also determinant in Chlorophyceae distribution. Guidelines are discussed for future, more robust and precise analyses of environmental algal DNA in mountain ecosystems and address green algae species distribution and dynamics in response to environmental changes.

Double-digest RAD-sequencing: do pre- and post-sequencing protocol parameters impact biological results?

Next-generation sequencing technologies have opened a new era of research in population genetics. Following these new sequencing opportunities, the use of restriction enzyme-based genotyping techniques, such as restriction site-associated DNA sequencing (RAD-seq) or double-digest RAD-sequencing (ddRAD-seq), has dramatically increased in the last decade. From DNA sampling to SNP calling, the laboratory and bioinformatic parameters of enzyme-based techniques have been investigated in the literature. However, the impact of those parameters on downstream analyses and biological results remains less documented. In this study, we investigated the effects of sevral pre- and post-sequencing settings on ddRAD-seq results for two biological systems: a complex of butterfly species (Coenonympha sp.) and several populations of common beech (Fagus sylvatica). Our results suggest that pre-sequencing parameters (i.e., DNA quantity, number of PCR cycles during library preparation) have a significant impact on the number of recovered reads and SNPs, on the number of unique alleles and on individual heterozygosity. In the same way, we found that post-sequencing settings (i.e., clustering and minimum coverage thresholds) influenced loci reconstruction (e.g., number of loci, mean coverage) and SNP calling (e.g., number of SNPs; heterozygosity) but had only a marginal impact on downstream analyses (e.g., measure of genetic differentiation, estimation of individual admixture, and demographic inferences). In addition, replication analyses confirmed the reproducibility of the ddRAD-seq procedure. Overall, this study assesses the degree of sensitivity of ddRAD-seq data to pre- and post-sequencing protocols, and illustrates its robustness when studying population genetics.

Speciation with gene flow: Evidence from a complex of alpine butterflies (Coenonympha, Satyridae).

Until complete reproductive isolation is achieved, the extent of differentiation between two diverging lineages is the result of a dynamic equilibrium between genetic isolation and mixing. This is especially true for hybrid taxa, for which the degree of isolation in regard to their parental species is decisive in their capacity to rise as a new and stable entity. In this work, we explored the past and current patterns of hybridization and divergence within a complex of closely related butterflies in the genus Coenonympha in which two alpine species, C. darwiniana and C. macromma, have been shown to result from hybridization between the also alpine C. gardetta and the lowland C. arcania. By testing alternative scenarios of divergence among species, we show that gene flow has been uninterrupted throughout the speciation process, although leading to different degrees of current genetic isolation between species in contact zones depending on the pair considered. Nonetheless, at broader geographic scale, analyses reveal a clear genetic differentiation between hybrid lineages and their parental species, pointing out to an advanced stage of the hybrid speciation process. Finally, the positive correlation observed between ecological divergence and genetic isolation among these butterflies suggests a potential role for ecological drivers during their speciation processes.

Unravelling the invasion history of the Asian tiger mosquito in Europe.

Multiple introductions are key features for the establishment and persistence of introduced species. However, little is known about the contribution of genetic admixture to the invasive potential of populations. To address this issue, we studied the recent invasion of the Asian tiger mosquito (Aedes albopictus) in Europe. Combining genome-wide single nucleotide polymorphisms and historical knowledge using an approximate Bayesian computation framework, we reconstruct the colonization routes and establish the demographic dynamics of invasion. The colonization of Europe involved at least three independent introductions in Albania, North Italy and Central Italy that subsequently acted as dispersal centres throughout Europe. We show that the topology of human transportation networks shaped demographic histories with North Italy and Central Italy being the main dispersal centres in Europe. Introduction modalities conditioned the levels of genetic diversity in invading populations, and genetically diverse and admixed populations promoted more secondary introductions and have spread farther than single-source invasions. This genomic study provides further crucial insights into a general understanding of the role of genetic diversity promoted by modern trade in driving biological invasions.

Body size determines soil community assembly in a tropical forest.

Tropical forests shelter an unparalleled biological diversity. The relative influence of environmental selection (i.e., abiotic conditions, biotic interactions) and stochastic-distance-dependent neutral processes (i.e., demography, dispersal) in shaping communities has been extensively studied for various organisms, but has rarely been explored across a large range of body sizes, in particular in soil environments. We built a detailed census of the whole soil biota in a 12-ha tropical forest plot using soil DNA metabarcoding. We show that the distribution of 19 taxonomic groups (ranging from microbes to mesofauna) is primarily stochastic, suggesting that neutral processes are prominent drivers of the assembly of these communities at this scale. We also identify aluminium, topography and plant species identity as weak, yet significant drivers of soil richness and community composition of bacteria, protists and to a lesser extent fungi. Finally, we show that body size, which determines the scale at which an organism perceives its environment, predicted the community assembly across taxonomic groups, with soil mesofauna assemblages being more stochastic than microbial ones. These results suggest that the relative contribution of neutral processes and environmental selection to community assembly directly depends on body size. Body size is hence an important determinant of community assembly rules at the scale of the ecological community in tropical soils and should be accounted for in spatial models of tropical soil food webs.

Diet shifts by adult flightless dung beetles Circellium bacchus, revealed using DNA metabarcoding, reflect complex life histories.

Life history changes may change resource use. Such shifts are not well understood in the dung beetles, despite recognised differences in larval and adult feeding ability. We use the flightless dung beetle Circellium bacchus to explore such shifts, identifying dung sources of adults using DNA metabarcoding, and comparing these with published accounts of larval dung sources. C. bacchus is traditionally considered to specialise on the dung of large herbivores for both larval and adult feeding. We successfully extracted mammal DNA from 151 adult C. bacchus fecal samples, representing 16 mammal species (ranging from elephants to small rodents), many of which are hitherto undescribed in the diet. Adult C. bacchus showed clear dung source preferences, especially for large herbivores inhabiting dense-cover vegetation. Our approach also confirmed the presence of cryptic taxa in the study area, and we propose that this may be used for biodiversity survey and monitoring purposes. Murid rodent feces were the most commonly fed-upon dung source (77.5%) for adult C. bacchus, differing markedly from the large and megaherbivore dung sources used for larval rearing. These findings support the hypothesis of life history-specific shifts in resource use in dung beetles, and reveal a hitherto unsuspected, but ecologically important, role of these dung beetles in consuming rodent feces. The differences in feeding abilities of the larval and adult life history stages have profound consequences for their resource use and foraging strategies, and hence the ecological role of dung beetles. This principle and its ecological consequences should be explored in other scarabaeids.

Differences in the fungal communities nursed by two genetic groups of the alpine cushion plant, Silene acaulis.

Foundation plants shape the composition of local biotic communities and abiotic environments, but the impact of a plant's intraspecific variations on these processes is poorly understood. We examined these links in the alpine cushion moss campion (Silene acaulis) on two neighboring mountain ranges in the French Alps. Genotyping of cushion plants revealed two genetic clusters matching known subspecies. The exscapa subspecies was found on both limestone and granite, while the longiscapa one was only found on limestone. Even on similar limestone bedrock, cushion soils from the two S. acaulis subspecies deeply differed in their impact on soil abiotic conditions. They further strikingly differed from each other and from the surrounding bare soils in fungal community composition. Plant genotype variations accounted for a large part of the fungal composition variability in cushion soils, even when considering geography or soil chemistry, and particularly for the dominant molecular operational taxonomic units (MOTUs). Both saprophytic and biotrophic fungal taxa were related to the MOTUs recurrently associated with a single plant genetic cluster. Moreover, the putative phytopathogens were abundant, and within the same genus (Cladosporium) or species (Pyrenopeziza brassicae), MOTUs showing specificity for each plant subspecies were found. Our study highlights the combined influences of bedrock and plant genotype on fungal recruitment into cushion soils and suggests the coexistence of two mechanisms, an indirect selection resulting from the colonization of an engineered soil by free-living saprobes and a direct selection resulting from direct plant-fungi interactions.

Genetic diversity and distribution differ between long-established and recently introduced populations in the invasive mosquito Aedes albopictus.

The Asian tiger mosquito Aedes albopictus, native to South-eastern Asia, is currently the most invasive mosquito in the world. The spatio-temporal dynamics of its expansion through the genetic characterization of invasive populations has been challenged so far by the limited number of genetic markers variable enough to infer the genetic structure in recently invaded areas. Here we applied the double-digest Restriction-site Associated DNA sequencing method (ddRADseq) to mosquitoes collected in two invaded areas, Reunion Island (12 localities) and Europe (18 localities). Analyses of genetic diversity, Bayesian clustering, Maximum Likelihood inference and isolation-by-distance tests based on 1561 genome-wide distributed Single Nucleotide Polymorphisms (SNPs) revealed that Reunion Island and Europe form two distinct genetic clusters, supporting no contemporary gene flow and suggesting two different and independent invasion histories. Long-established populations (Reunion Island) were more genetically diverse than recently introduced European populations. The largest part of genetic variance was found at the intra-individual level (>85%) and most FIS values were positive, suggesting inbreeding at the local scale. The two invaded areas showed contrasting patterns of genetic structure. Significant isolation-by-distance was found among Reunion Island populations, suggesting that these populations are at the drift-migration equilibrium. In contrast, long-distance human-assisted transport is probably the main dispersal mechanism in Europe.

At the Origin of a Worldwide Invasion: Unraveling the Genetic Makeup of the Caribbean Bridgehead Populations of the Dengue Vector Aedes aegypti.

Human-driven global environmental changes have considerably increased the risk of biological invasions, especially the spread of human parasites and their vectors. Among exotic species that have major impacts on public health, the dengue fever mosquito Aedes aegypti originating from Africa has spread worldwide during the last three centuries. Although considerable progress has been recently made in understanding the history of this invasion, the respective roles of human and abiotic factors in shaping patterns of genetic diversity remain largely unexplored. Using a genome-wide sample of genetic variants (3,530 ddRAD SNPs), we analyzed the genetic structure of Ae. aegypti populations in the Caribbean, the first introduced territories in the Americas. Fourteen populations were sampled in Guyane and in four islands of the Antilles that differ in climatic conditions, intensity of urbanization, and vector control history. The genetic diversity in the Caribbean was low (He = 0.14-0.17), as compared with a single African collection from Benin (He = 0.26) and site-frequency spectrum analysis detected an ancient bottleneck dating back ∼300 years ago, supporting a founder event during the introduction of Ae. aegypti. Evidence for a more recent bottleneck may be related to the eradication program undertaken on the American continent in the 1950s. Among 12 loci detected as FST-outliers, two were located in candidate genes for insecticide resistance (cytochrome P450 and voltage-gated sodium channel). Genome-environment association tests identified additional loci associated with human density and/or deltamethrin resistance. Our results highlight the high impact of human pressures on the demographic history and genetic variation of Ae. aegypti Caribbean populations.

Is There Any Evidence for Rapid, Genetically-Based, Climatic Niche Expansion in the Invasive Common Ragweed?

Climatic niche shifts have been documented in a number of invasive species by comparing the native and adventive climatic ranges in which they occur. However, these shifts likely represent changes in the realized climatic niches of invasive species, and may not necessarily be driven by genetic changes in climatic affinities. Until now the role of rapid niche evolution in the spread of invasive species remains a challenging issue with conflicting results. Here, we document a likely genetically-based climatic niche expansion of an annual plant invader, the common ragweed (Ambrosia artemisiifolia L.), a highly allergenic invasive species causing substantial public health issues. To do so, we looked for recent evolutionary change at the upward migration front of its adventive range in the French Alps. Based on species climatic niche models estimated at both global and regional scales we stratified our sampling design to adequately capture the species niche, and localized populations suspected of niche expansion. Using a combination of species niche modeling, landscape genetics models and common garden measurements, we then related the species genetic structure and its phenotypic architecture across the climatic niche. Our results strongly suggest that the common ragweed is rapidly adapting to local climatic conditions at its invasion front and that it currently expands its niche toward colder and formerly unsuitable climates in the French Alps (i.e. in sites where niche models would not predict its occurrence). Such results, showing that species climatic niches can evolve on very short time scales, have important implications for predictive models of biological invasions that do not account for evolutionary processes.

Hybridization promotes speciation in Coenonympha butterflies.

Hybridization has become a central element in theories of animal evolution during the last decade. New methods in population genomics and statistical model testing now allow the disentangling of the complexity that hybridization brings into key evolutionary processes such as local adaptation, colonization of new environments, species diversification and extinction. We evaluated the consequences of hybridization in a complex of three alpine butterflies in the genus Coenonympha, by combining morphological, genetic and ecological analyses. A series of approximate Bayesian computation procedures based on a large SNP data set strongly suggest that the Darwin's Heath (Coenonympha darwiniana) originated through hybridization between the Pearly Heath (Coenonympha arcania) and the Alpine Heath (Coenonympha gardetta) with different parental contributions. As a result of hybridization, the Darwin's Heath presents an intermediate morphology between the parental species, while its climatic niche seems more similar to the Alpine Heath. Our results also reveal a substantial genetic and morphologic differentiation between the two geographically disjoint Darwin's Heath lineages leading us to propose the splitting of this taxon into two different species.

Upscaling the niche variation hypothesis from the intra- to the inter-specific level.

The "niche variation hypothesis" (NVH) predicts that populations with wider niches should display higher among-individual variability. This prediction originally stated at the intra-specific level may be extended to the inter-specific level: individuals of generalist species may differ to a greater extent than individuals of a specialist species. We tested the NVH at intra- and inter-specific levels based on a large diet database of three large herbivore feces collected in the field and analyzed using DNA metabarcoding. The three herbivores (roe deer Capreolus capreolus, chamois Rupicapra rupicapra and mouflon Ovis musimon) are highly contrasted in terms of sociality (solitary to highly gregarious) and diet. The NVH at the intraspecific level was tested by relating, for the same population, diet breadth and inter-individual variation across the four seasons. Compared to null models, our data supported the NVH both at the intra- and inter-specific levels. Inter-individual variation of the diet of solitary species was not larger than in social species, although social individuals feed together and could therefore have more similar diets. Hence, the NVH better explained diet breadth than other factors such as sociality. The expansion of the population niche of the three species was driven by resource availability, and achieved by an increase in inter-individual variation, and the level of inter-individual variability was larger in the generalist species (mouflon) than in the specialist one (roe deer). This mechanism at the base of the NVH appears at play at different levels of biological organization, from populations to communities.

Fifty thousand years of Arctic vegetation and megafaunal diet.

Although it is generally agreed that the Arctic flora is among the youngest and least diverse on Earth, the processes that shaped it are poorly understood. Here we present 50 thousand years (kyr) of Arctic vegetation history, derived from the first large-scale ancient DNA metabarcoding study of circumpolar plant diversity. For this interval we also explore nematode diversity as a proxy for modelling vegetation cover and soil quality, and diets of herbivorous megafaunal mammals, many of which became extinct around 10 kyr bp (before present). For much of the period investigated, Arctic vegetation consisted of dry steppe-tundra dominated by forbs (non-graminoid herbaceous vascular plants). During the Last Glacial Maximum (25-15 kyr bp), diversity declined markedly, although forbs remained dominant. Much changed after 10 kyr bp, with the appearance of moist tundra dominated by woody plants and graminoids. Our analyses indicate that both graminoids and forbs would have featured in megafaunal diets. As such, our findings question the predominance of a Late Quaternary graminoid-dominated Arctic mammoth steppe.

Multi-stemmed trees of Nothofagus pumilio second-growth forest in Patagonia are formed by highly related individuals.

BACKGROUND AND AIMS: Multi-stemmed trees (tree clusters) in Nothofagus pumilio, a dominant tree species in Patagonia, are very uncommon and are restricted to the edge of second-growth forests following human-provoked fires. No vegetative reproduction has been reported so far. The genetic structure of multi-stemmed trees of this species was investigated and it was hypothesized that genets within a cluster were more closely related than average in the population. METHODS: Fifteen clusters (composed of at least three purported stems) and 15 single trees were sampled at the edge of a second-growth forest and genotyped using two amplified fragment length polymorphism (AFLP) primer pairs. We obtained 119 polymorphic markers that allowed clonality to be determined, together with sibship structure and relatedness among samples. KEY RESULTS: Clonality was detected in seven clusters but all clusters had at least two different genotypes. Full sibs were found exclusively within clusters and in all clusters. Within a cluster, stems that were not identified as full sibs were often half sibs. Relatedness values for the full sibs and half sibs were higher than the theoretical values of 0·5 and 0·25 but the relatedness between clusters was very low. CONCLUSIONS: Tree clusters that are merged at the edge of the second-growth forest of N. pumilio are composed of stems of the same genotype and of other genotypes that are highly related (but not always). It is suggested that this peculiar genetic structure results from a combination of several causes, including selection for merging of related individuals.

Soil sampling and isolation of extracellular DNA from large amount of starting material suitable for metabarcoding studies.

DNA metabarcoding refers to the DNA-based identification of multiple species from a single complex and degraded environmental sample. We developed new sampling and extraction protocols suitable for DNA metabarcoding analyses targeting soil extracellular DNA. The proposed sampling protocol has been designed to reduce, as much as possible, the influence of local heterogeneity by processing a large amount of soil resulting from the mixing of many different cores. The DNA extraction is based on the use of saturated phosphate buffer. The sampling and extraction protocols were validated first by analysing plant DNA from a set of 12 plots corresponding to four plant communities in alpine meadows, and, second, by conducting pilot experiments on fungi and earthworms. The results of the validation experiments clearly demonstrated that sound biological information can be retrieved when following these sampling and extraction procedures. Such a protocol can be implemented at any time of the year without any preliminary knowledge of specific types of organisms during the sampling. It offers the opportunity to analyse all groups of organisms using a single sampling/extraction procedure and opens the possibility to fully standardize biodiversity surveys.

Selection criteria for scoring amplified fragment length polymorphisms (AFLPs) positively affect the reliability of population genetic parameter estimates.

A reliable data set is a fundamental prerequisite for consistent results and conclusions in population genetic studies. However, marker scoring of genetic fingerprints such as amplified fragment length polymorphisms (AFLPs) is a highly subjective procedure, inducing inconsistencies owing to personal or laboratory-specific criteria. We applied two alternative marker selection algorithms, the newly developed script scanAFLP and the recently published AFLPScore, to a large AFLP genome scan to test how population genetic parameters and error rates were affected. These results were confronted with replicated random selections of marker subsets. We show that the newly developed marker selection criteria reduced the mismatch error rate and had a notable influence on estimates of genetic diversity and differentiation. Both effects are likely to influence biological inference. For example, genetic diversity (HS) was 29% lower while genetic differentiation (FST) was 8% higher when applying scanAFLP compared with AFLPScore. Likewise, random selections of markers resulted in substantial deviations of population genetic parameters compared with the data sets including specific selection criteria. These randomly selected marker sets showed surprisingly low variance among replicates. We conclude that stringent marker selection and phenotype calling reduces noise in the data set while retaining patterns of population genetic structure.

Tracking genes of ecological relevance using a genome scan in two independent regional population samples of Arabis alpina.

Understanding the genetic basis of adaptation in response to environmental variation is fundamental as adaptation plays a key role in the extension of ecological niches to marginal habitats and in ecological speciation. Based on the assumption that some genomic markers are correlated to environmental variables, we aimed to detect loci of ecological relevance in the alpine plant Arabis alpina L. sampled in two regions, the French (99 locations) and the Swiss (109 locations) Alps. We used an unusually large genome scan [825 amplified fragment length polymorphism loci (AFLPs)] and four environmental variables related to temperature, precipitation and topography. We detected linkage disequilibrium among only 3.5% of the considered AFLP loci. A population structure analysis identified no admixture in the study regions, and the French and Swiss Alps were differentiated and therefore could be considered as two independent regions. We applied generalized estimating equations (GEE) to detect ecologically relevant loci separately in the French and Swiss Alps. We identified 78 loci of ecological relevance (9%), which were mainly related to mean annual minimum temperature. Only four of these loci were common across the French and Swiss Alps. Finally, we discuss that the genomic characterization of these ecologically relevant loci, as identified in this study, opens up new perspectives for studying functional ecology in A. alpina, its relatives and other alpine plant species.

Evolution and taxonomy of the wild species of the genus Ovis (Mammalia, Artiodactyla, Bovidae).

New insights for the systematic and evolution of the wild sheep are provided by molecular phylogenies inferred from Maximum parsimony, Bayesian, Maximum likelihood, and Neighbor-Joining methods. The phylogeny of the wild sheep was based on cytochrome b sequences of 290 samples representative of most of the sub-species described in the genus Ovis. The result was confirmed by a combined tree based on cytochrome b and nuclear sequences for 79 Ovis samples representative of the robust clades established with mitochondrial data. Urial and mouflon, which are either considered as a single or two separate species, form two monophyletic groups (O. orientalis and O. vignei). Their hybrids appear in one or the other group, independently from their geographic origin. The European mouflon O. musimon is clearly in the O. orientalis clade. The others species, O. dalli, O. canadensis, O. nivicola, and O. ammon are monophyletic. The results support an Asiatic origin of the genus Ovis, followed by a migration to North America through North-Eastern Asia and the Bering Strait and a diversification of the genus in Eurasia less than 3 million years ago. Our results show that the evolution of the genus Ovis is a striking example of successive speciation events occurring along the migration routes propagating from the ancestral area.